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Abstract 

A calculation of the blackbody radiation shift of the clock transition is performed. The po- 
larizabilities of the 2s^ ^S"^, 2s2p ^P°, and 2s2p "^P" states are computed using the configuration 
interaction method with an underlying semi-empirical core potential. The recommended dipole po- 
^ \ larizabilities are 9.64(3) Cg, 7.78(3) aj] and 16.55(5) Og respectively. The derived frequency shift for 

^ the 2^2 ^S^^ 2s2p ^Pg transition at 300 K is 0.0160(5) Hz. The dipole polarizabilities agree with 

an earlier relativistic calculation (Safronova et al. Phys. Rev. Lett. 107 143006 (2011)) to better 
than 0.2%. Quadrupole and octupole polarizabilities and non-adiabatic multipole polarizabilities 
are also reported. 
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I. INTRODUCTION 



Recent advances in laser control of atoms and ions have lead to major improvements in 
the precision of optical frequency standards [l|, [2|. These improvements are expected to 
result in a new definition of the second [3| . Indeed, an optical clock using the A1+ ion using 
quantum logic technology has been developed with a fractional frequency uncertainty of 
8.6 X 10~^® 1^. This uncertainty is equivalent to a drift of 1 second in 3.7 x 10^ years. 

The ultra-high precision achieved by these optical frequency standards means they are 
sensitive to very small environmental influences. One of these influences is the blackbody 
radiation (BBR) emitted by the apparatus in which the atomic or ionic clock is enclosed. 
This BBR radiation, by means of the AC Stark effect changes the energies of the two states 
of the clock transition, and this can alter the frequency of the atomic clock js], [6|. 

This BBR-shift is in principle one of the largest sources of systematic error in these clocks 
The BBR-shift (in Hz) can be written 



Az/BBR = 6.579684 x lO^^ (AE„ - A^i^^er) , (1) 

where the electric dipole (El) induced BBR energy shift of an atomic state can be approxi- 
mately calculated as [l2| 

Ai?^-l(a7r)V(0)T^ (2) 

The dipole polarizability of the relevant quantum state is ai and T is the temperature. 
Knowledge of the dipole polarizabilities permits a temperature dependent BBR correction 
to be made to the clock. The uncertainty in the El BBR shift can be written 

«(A.eBK) = A.BBK + . (3) 

Calculations of the B"*" (2s^ ^Sq - 2s2p '^Pq) clock transition have previously been made 



13( 1 using a relativistic configuration interaction (CI) calculation to account for valence 
correlations while an all-order many-body perturbation theory approach is used to account 
for core and core-valence correlations. The paper reported the dipole polarizabilities and 
demonstrated that the clock transition for this ion had a relatively small 300 K BBR shift 
of 0.0159 Hz. The present manuscript confirms this result and extends the dataset for B"*" 
to encompass higher order polarizabilities. Calculations are performed using the CI method 
with a semi-empirical core-polarization potential to encompass core-valence correlations. 

II. METHODOLOGY 

The CI calculations used to generate the physical and pseudo states were similar in 
style to those used previously to determine the dispersion parameters and polarizabilities of 



a number of two electron systems |14|-|17|| . The Hamiltonian for the two active electrons is 
written 

^ = 11 (-^^' + ^dir(r.) + V;xc(r.) + l^pi(r,)^ 

+ ■^"p2(ri,r2) + — . (4) 
r 12 
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The direct, Vdir, and exchange, Vexc, interactions of the valence electrons with the Hartree- 
Fock (HF) core were calculated exactly. The Is^ core wave function was taken from a HF 
calculation of the B^"*" ground state using a Slater type orbital (STO) basis. The ^-dependent 
polarization potential, Vpi, was semi-empirical in nature with the functional form 



^Pi(r) = -E^'^?|^IM(M- (5) 



Cm 



2r^ 



The coefficient, ttcore, is the static dipole polarizabihty of the core and gj{r) = 1 — exp(- 
/ pf) is a cutoff function designed to make the polarization potential finite at the origin. 
The cutoff parameters, p^, were tuned to reproduce the binding energies of the B^"*" ns 
ground state and the np, n d and nf excited states. The core polarizabihty was chosen to 
be acore = 0.019644 (isl, [T^- The cutoff parameters for £ = ^ 3 were 0.6835, 0.6899, 



0.8874 and 2.945 ao respectively. 

To get more accurate energy levels and polarizabilities, it is essential to include a two 
body polarization term, Vp2-, in the Hamiltonian. The polarization of the core by one electron 
is influenced by the presence of the second valence electron. Omission of the two-body term 
would typically result in a 2s^ state that would be too tightly bound. A discussion of the 
importance of the two body polarization potential can be found in 23|. The two body 



polarization potential is adopted in the present calculation with the form 

Vp2{ju Tj) = — ^(ri ■ rj)^p2(ri)^p2(rj) , (6) 

where gp2 has the same functional form as gi{r). The cutoff parameter for gp2{i^) was chosen 
as 0.6867 ao, the average of po and pi (the p2 and ps cutoff parameters are influenced by 
finite nuclear mass effects, and thus they were not used in determining the cutoff parameter 
for Vp2)- Use of 0.6867 Oq for the two-body cutoff parameter resulted in energies that were 
close to the experimental binding energies for most of the lowest lying states of B+. Some 
small adjustments to the pi, described later, were made later to further improve agreement 
with the experimental B^ spectrum. The approach to solve the Schrodinger equation is 
termed as conflguration interaction plus core polarization (CICP). 

There were a total of 163 valence orbitals with a maximum orbital angular momentum 
oi £ = 5. The radial dependence of the orbitals were described by a mixture of STOs and 
Laguerre type orbitals (LTOs) [13]. The number of active orbitals for £ = — ?■ 5 were 32, 
32, 30, 25, 25, and 19 respectively. Some i = valence orbitals were generated from the 
STOs used for the core. All the other orbitals were written as LTOs due to their superior 
linear dependence properties when compared with STO basis sets. The use of the large 
orbital basis resulted in wave functions and energies for the low-lying states that were close 
to convergence. 

The length of the CI expansions for the different states of B"*" ranged from 2000-5000. 
Some small changes were made to the pe values that were originally tuned to the B^"*" 
spectrum to improve the agreement of the B"*" energies with experiment. The oscillator 



strengths were computed with operators that included polarization corrections |14j . l2ll . |22 



The cutoff parameter in the polarization correction to dipole operator was 0.6867 Oq. 
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TABLE I: Theoretical and experimental energy levels (in Hartree) for some of the low-lying 
states of the B+ and ions. The energies are given relative to the energy of the B'^+ core. The 
experimental energies for the multiplet states are averages with the usual (2J+1) weighting factors. 
The experimental data were taken from the National Institute of Standards and Technology [23^]. 



State Present Experiment 





B2+ 




2s ^S*" 


-1 393924 


— 1 393924 


2p 


-1 173483 


— 1 173483 


3s ^S*" 


-0 572792 


-0 572863 


3p 


-0.514642 


-0 514743 


3d '^D^ 


-0 500561 


-0 500561 

Vy a w vy vy w -L 


4s ^S*" 


-0 310856 


—0 310891 


4o ^P° 


-0 287444 


—0 287498 




-0 281527 


-0.281529 


4f 


-0.281269 


-0.281269 




B+ 




2s2 iS^ 


—2 318347 


—2 318347 


ZsZp r 






2s2p ip" 


-1.983927 


-1.983927 


2p2 3pe 


-1.867605 


-1.867605 


2p2 


-1.851947 


-1.851947 


2p2 is^ 


-1.736606 


-1.736679 


2s3s ^S'^ 


-1.727053 


-1.727053 


2s3s ^S^ 


-1.691092 


-1.691293 


2s3p ^P" 


-1.662237 


-1.662269 


2s3j9 ip° 


-1.661828 


-1.661765 


2s3d ^D^ 


-1.631961 


-1.631961 


2s3d ^D^ 


-1.613484 


-1.613545 



III. RESULTS AND DISCUSSION 
A. Energy levels 

The energy levels of the present calculations are given in Table |T] and compared with 
experiment. The biggest discrepancy for the ion was 10""^ a.u.. The cut-off parameters 
of the the polarization potential were tuned to reproduce the experimental binding energies 
of the lowest states of each symmetry. 

Small adjustments to the cut-off parameters were made for the calculations of the 
states. For example, the value of po was reset to 0.7064 ao for the calculation of the states of 
the ^S*^ symmetry. The value of po was fixed by requiring that the theoretical and experimen- 
tal energies for the state be the same. Other fine tunings of the cut-off parameters were 
made for all symmetries. The biggest discrepancy between theoretical and experimental 
energies occurs for the ^S^ symmetry and is only 2 x 10^^ a.u.. 
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The agreement between the theoretical and experimental energy levels is sufficiently close 
to discount the possibility that energy level considerations might make a significant contri- 
bution to the uncertainty in the radial matrix elements. 



B. Oscillator strengths of low-lying transitions 



TABLE II: Absorption oscillator strengths for various dipole transition lines of the B"*" and B-^+ 
ions. The experimental energy differences were used in the calculation of the CICP oscillator 
strengths. 



Transition 


CICP 


BCICP 


MCHF-BP 


MCHF 


CI 


Other Theory 


Experiment 












B2+ 










2s-*-2p 


0.36360 






0.36370 [24| 




0.36389 [25] 


0.363243 


' [261 


0.35(2) [27] 


2s->3p 


0.15333 






0.15346 [24| 




0.15376 [28] 






0.15(1) [27] 


2s— >4p 


0.04969 










0.04981 [28] 








2p—^3s 


0.04640 






0.04636 124] 










0.05(1) [27] 


2p^3d 


0.63801 






0.63803 [24] 


B+ 








0.62(6) [27] 


2s2 lS'=->2s2p lp° 


0.99907 


1.002 


[29] 


1.001 [30| 


0.9976(22) [31] 


0.9997 [32] 


1.0012' 


[33] 


0.98(8) [34] 












0.999(5) [35] 


1.005 [36] 


1.0012'= 


[37] 


0.71(5) [27] 














0.9998'' [38] 






0.98(6) [39] 














1.0028= [13] 








2s2 lS'=^2s3p lp° 


0.10959 


0.108 


[29] 


0.1087 [30] 


0.1093(3) [40] 










2s2p lp°^2p2 iD^ 


0.16195 


0.162 


[29] 


0.1621 [30] 


0.1608(44) [31] 


0.1625 [32] 






0.192(9) [41] 




















0.114(6) [42] 


2s2p ip°^2s3d 


0.51545 


0.514 


[29] 


0.5161 [30J 




0.5199 [32] 






0.49(2) [27] 


2s2p ip°^2p2 iS'= 


0.22591 


0.227 


[29] 


0.2259 [30] 


0.2257(38) [31] 


0.2264 [32] 






0.24(2) [34] 




















0.20(1) [27] 




















0.163(11) [41] 


2s2p lp°^2s3s ^3= 


0.00008 






0.00019 [30] 




0.00007 [32] 






0.039(2) [27] 


2s2p ^P°^2p2 3pe 


0.34298 


0.365 


[43] 


0.34292 [^] 


0.3427(2) [31] 


0.3427 [32] 






0.34(3) [34] 




















0.32(2) [27] 


2s2p ^P°^2s3s ^'S'^ 


0.06377 






0.06401 [30| 












2s2p ^P°^2s3d ^D" 


0.47627 


0.473 


[29] 


0.47597 [30] 










0.49(2) [27] 


2p2p 3p<!^2p3d 3d° 


0.62300 


0.310 


[43] 















"Hylleraas— type variational metliod. 
^Multi-Configuration Dirac-Fock (MCDF) method. 
"^Relativistic MBPT calculation. 
'^Relativistic CI calculation with MBPT theory. 

'^Relativistic CI calculation with all order MBPT theory. Calculated with theoretical energy differences. 



The oscillator strengths for the transitions between the low manifolds states are listed in 
Tables HTl The absorption oscillator strength from state '^i to state '^j is calculated according 



to the identity [IJ, |44 



(fe) 



(2fc + l)(2L, + l) 



(7) 



In this expression, Sji = {Ej — Ei) is the energy difference between the initial state and 
final state, while k is the multipolarity of the transition, and C'^(f) is a spherical tensor. 
Experimental energy differences were used for the calculation of oscillator strengths. 
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There have been very many calculations performed of the ener gy l evels and oscillator 



33, 35 



40, 43, 47-51 . Not all of the 



strengths for 3^+ l28|, |45l |46 



and B+ [29- 

theoretical calculations were tabulated. Table HTl gives the reported results of the calculations 
that are deemed to be the most accurate or of particular relevance to present calculations. 

For ion /-values are given the Multi- Configuration Hartree-Fock calculation with 
Briet-Pauli corrections (MCHF-BP) The present calculations agree with the MCHF-BP 
values to an accuracy of 0.0001. While the present calculations are ostensibly nonrelativistic, 
they implicitly include relativistic corrections since the energies are tuned to experimental 
values. The 2s — 2p oscillator strength computed with the Hylleraas method [2^ is close to 
the non-relativist ic limit, but the Hylleraas calculation omits any relativistic effects and the 
Hylleraas energy difference for the 2s — )■ 2p transitions is 0.22016 a.u. which is about 0.1% 
smaller than the experimental energy difference. The full core plus correlation calculation 
listed in the CI column is a variant of the configuration interaction approach. 



[251, |28|, 143 



There is one previous calculation for B+ that is very similar in concept to the present 



ii. The 



methodology. That was a CI calculation with a semi-empirical core potential [291 . 
major distinction was the adoption of a B-spline basis so this calculation is abbreviated as 
BCICP in Table [Tll With a few exceptions, the CICP and BCICP oscillator strengths agree 
to about 1%. When the BCICP oscillator strengths are different from the present values, one 
also finds the BCICP oscillator strengths also disagreeing with the MCHF-BP B^+ oscillator 
strengths 

There is also better than 0.3% agreement of the CICP calculation with MCHF oscillator 
strengths with two exceptions. The MCHF oscillator strength 31| for the 2s2p — )• 2p^ 
^D*^ transition is about 1% smaller than the CICP oscillator strength. The MCHF oscillator 
strength however is about 1% smaller that the BCICP and MCHF-BP oscillator strengths. 
There is also agreement at better than 1% level with a CI calculation [3^ except for the 
case of the 2s2p ^P° — )■ 2s3s ^S^ transition which has a very small oscillator strength. 

There have been two calculations which combine relativistic CI calculations with many 



body perturbation theory (CI+MBPT) to represent the core-valence interaction [13|, |38 
These only gave the oscillator strength for the 2s^ ^S*^ to 2s2p ^P° transition. The total 
range between the CICP oscillator strengths and two CI+MBPT oscillator strengths is less 
than 0.4%. The agreement of the CICP oscillator strengths with another two relativistic 

calculations that are the MCDF calculation of ^ ^ 

0.2% level. " 

Some experimental oscillator strength measurements [2 
Table [TTl for completeness. The precision of the experimental data is not as high as many of 
the theoretical oscillator strengths. 



33 and MBPT calculation of 37 is also at 



H, m, SH Q are also listed in 



C. Scalar and tensor polarizabilities 

This analysis is done under the premise that spin-orbit effects are small and the radial 
parts of the wave functions are the same for the states with different J. 

All the polarization parameters reported here are calculated using their respective os- 
cillator strength sum rules. The multipole oscillator strengths f-f^ are defined in Eq. ([7]). 



Then the adiabatic multipole polarizabilities from the state i are written as [52 

J 1 
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TABLE III: The pseudo-oscillator strength distribution for the core B^+. The energy shift pa- 
rameter ej and the adiabatic (ocore) and non-adiabatic (/3core) core polarizabilities from Hylleraas 
calculations 



18( 1 are also displayed. The numbers in the square brackets denote powers of 10. 



dipole -16.67592 1.0 1.9644[-2] 1.1243[-3] 

-7.89382 1.0 

quadrupole -21.91592 0.28537 3.4266[-3] 1.5237[-4] 

-10.14212 0.28537 

octupole -22.51592 0.15844 1.5216[-3] 5.9751[-5] 

-11.44722 0.15844 



Related sum rules such as the non-adiabatic multipole polarizability and Sfc(-4) are given 
as [U] 

J !■ 



1 



3 



and l5i 

M-4) = E7F- (10) 



The 5'fc(— 4) sum rule gives the lowest order frequency dependent component to the dynamic 
polarizability through the relation 

= «fe(0) + a;2^fc(-4) + .... (11) 

States with a non-zero angular momentum will also have a tensor polarizability js], |53| . 
For a state with angular momentum Lq{Jq), this is defined as the polarizability of the 
magnetic sub- level with M = Lq{M = Jq). The total polarizability is written in terms of 
both a scalar and tensor polarizability. The scalar polarizability represents the average shift 
of the different M levels while the tensor polarizability gives the differential shift. 

This tensor polarizability can be expressed in terms of /-value sum rules. For an Lq = 1 
initial state, one can write the tensor polarizability for a dipole field as 

/ fon 1 fon , 1 fon\ 

V „ ,,L„=0 ^"0 ^ n,L„=l ^"0 n,L„=2 ^ 



If the initial state is a Lq = 2 state, one can use the expressions in [53| and get the /-value 
sum 

/ fon fon , 2 fon\ „x 

a2,LoLo = -( 2^ 2^ TT+n 2^ Tt] ■ (13) 

^ n,L„=l "0 n,L„=2 "0 n,L„=3 ^ 

The core does not make a contribution to the tensor polarizability since it has an equal 
impact on all the different M-levels. 

The development above is for LS coupled states, but it is common to give the tensor 
polarizability for LSJ states. These can be related to the LS states by geometric factors 
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TABLE IV: The polarizabilities of some low lying states of B^^ ion. The scalar adiabatic polar- 
izabilities (a^) are listed along with some non-adiabatic (/3fc) and tensor ict'2LQLo ) polarizabilities. 
All the polarizabilities are calculated using the experimental energies. The dipole polarizabilities 
from accurate CI calculations 25j, |5J] are displayed for comparison. The polarizabilities are in 
atomic units. 



State 


ai 


/3i 


2,LoLo 


a2 


/32 


2s 


7.8460 
7.847 [25] 
7.85 [54] 


17.137 




7.0963 


3.8719 


2p 2p° 


-0.56938 


6.9896 


2.1659 


5.6105 


3.1396 


3s 2S'= 


182.94 


1558.1 




1539.0 


9096.2 


3p 2p« 


312.04 


13153 


20.605 


643.82 


1374.8 


3d 20*= 


-191.26 


7616.0 


208.16 


10.437 


2170.3 



as 



30.181 



48.761 

14598 
65846 
-24466 



arising from the application of Racah algebra. The scalar polarizabilities for the different J 
levels are the same (if spin-orbit splitting is neglected) and equal to the scalar polarizability 
in the L representation. The tensor polarizabilities between the L and J representations 
can be related using the expressions of [s^]. When Lq = 1 and Jo = one finds a2,JoJo = 
while Jo = 1 case gives a2,JoJo = ^ "2,LoLo- 



1. Core polarizabilities 

The energy distribution of the oscillator strengths originating from core excitations was 
estimated using a semi-empirical technique [13] • This approach utilizes /-value sum rules 
to construct the pseudo-oscillator strength distributions, 

<^k,core = |~T^ ; (14) 

iGcore ^ 

where Ni is the number of electrons in a core orbital, is an energy shift parameter. The 
energy shift parameter was chosen so that Eq. f ll4p reproduces accurate estimates of the 
adiabatic and non-adiabatic core polarizabilities determined by close to exact calculations 
for dipole, quadrupole and octupole transitions [l8[. 

The present calculated pseudo-oscillator strength distributions are given in Table lllli 
They can be used in the determination of the dynamic polarizabilities and the long range 
van der Waals coefficients of the and B"*" ions with other atoms. 



2. The 5+ and B^^ polarizabilities 

Tables IIVI and |V] give the scalar adiabatic multipole polarizabilities of the lowest five 
states of the B^"*" ion and the lowest three states of the B+ ion. The tensor polarizabilities 
and non-adiabatic polarizabilities as well as the related sum rules 5'fc(— 4) of some states are 
also listed. The energies of the lowest lying states (i.e. those in Table [T]) were adjusted to 
be the same as the experimental energies for the polarizability calculations. 
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TABLE V: The polarizabilities of the lowest three states of the B+ ion. The scalar adiabatic 
polarizabilities (a^) are listed along with some non-adiabatic (/3fc) and tensor {ct'^l^gig ) polariz- 
abilities. Values for sum rules, 5i(— 4), are also presented. All these values are calculated using 
the experimental energies. The dipole polarizabilities from the relativistic calculation of H are 
displayed for comparison. All polarizabilities are in atomic units. 



State 


ai 


Si{-4) 


/3i 


2,LoL() 


a2 


/32 


"3 




9.6442 


80.891 


13.757 




27.138 


20.631 


147.01 




9.624 Cl+aU [13] 














2s2p 3p° 


7.7798 


66.320 


10.737 


1.4613 


25.011 


17.940 


220.98 




7.772 Cl+aU [13j 














2s2p ^P" 


16.554 


603.57 


54.587 


-2.1960 


44.757 


44.266 


541.44 



The present CICP dipole polarizability for the B^^ ground state is slightly smaller than 
the polarizability of two very accurate CI type calculations [25|, HI. The differences do not 
exceed 0.004 Oq. The CI calculations are non- relativistic and are expected to be slightly 
larger than the actual polarizability [HH']. A comparison for the iso-electronic ion Be"^ can be 
used to estimate an uncertainty in the B^"*" 2s states dipole polarizability. A previous CICP 
calculation gave a dipole polarizability of 24.493 Oq which is very close to the recommended 
value of 24.489(4) Oq jssf. Assigning an uncertainty of 0.1% to the B^^ ground state polar- 
izability would seem to be justified. The uncertainties in the polarizabilities for the excited 
states are expected to be of the same order as that of the ground state except for the case 
of the 2p state where considerable cancellations occur in the oscillator strength sum rule. 

Table |V] gives the polarizabilities of the B^ states. The only other calculation of the 
polarizabilities for these states is a recent CI-I-MBPT calculation [IS]. The CI-I-MBPT 
calculation gave a polarizability for the 2s'^ ground state that is 0.2% smaller than the 
present CICP calculation. The difference for the 2s2p state is 0.1%. 

A rough estimate of the uncertainties in the B"*" polarizabilities is possible by reference 



to similar calculations for the Si^"*" ground state 16|. A CICP calculation gave 11.688 Oq, a 



revised analysis of a resonant excitation Stark ionization spectroscopy (RESIS) experiment 
gave 11.669 al [fliej , and a CI+MBPT calculation gave 11.670(13) al (57|. The comparison 



between the RESIS and the CICP polarizabilities suggested that a conservative estimate of 
the ground state Ss^ ^S"^ polarizability was 0.25% while that for the excited states was 0.5%. 



Comparisons for A1+ between CICP polarizabilities [17[ and CI-I-MBPT calculations [13 
reveal differences between the two calculations that do not exceed 0.4%. The CI+MBPT 
calculation uses theoretical differences in the calculation of the A1+ polarizabilities and 
overestimates the 3^^ ^S*^ - 3s3p '^Pq energy difference by 0.14%. The replacement of the 
theoretical energy differences by the experimental energy differences would reduce the dif- 
ference between the CICP and CI+MBPT calculation to less than 0.3%. The analysis for 
the Al"*" system suggests that an uncertainty of 0.3% should be assigned to the polarizability 
of the 2^2 ^S^ state. So the final recommended 2^^ ^S^ dipole polarizability is 9.64(3) Oq. 
Assuming the 2s2p ^Pq state has the same uncertainty, the final CICP dipole polarizability 
is 7.78(3) Oq. The scalar dipole polarizability for 2s2p ^P° state was 16.55(5) assuming 
the same relative uncertainty. 

Table |VT] gives a breakdown of the different contributions to the 25^ ^S*^, 2s2p ^Pq, and 
2s2p ^P° dipole polarizabilities. About 90% of the 2s'^ ^S'^ polarizability comes from the 
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TABLE VI: Breakdown of the contributions to the dipole polarizabihties of the B+ clock transition 
states. The 5ai column gives the contribution from the indicated transition class. The ^ ai column 



gives the accumulated sum. 


The final polarizabihties are given in bold-face. 




Transition(s) 


5ai 








2s2 iS*^ state 






2s2p 


8.9333 


8.9333 




2s3p ^P° 


0.2542 


9.1875 


2s2 ^ 


nP 


0.4370 


9.6245 


Core 




0.01964 


9.6441 






2s2p 3po state 




2s2p 


2s3s ^S*" 


0.3596 


0.3596 


2s2p ^P" 


n5 ^S*^ 


0.1093 


0.4689 


2s2p 


-)• 2p2 ^P"^ 


4.3573 


4.8262 


2s2p ^P° 


nP ^P"^ 


0.0750 


4.9012 


2s2p ^P° 


2s3d ^D^ 


1.7873 


6.6885 


2s2p ^P° 




1.0717 


7.7602 


Core 




0.01964 


7.7798 






2s2j5 ip° state 




2s2p ^P° 


2s2 


-2.9778 


-2.9778 


2s2p 


2p2 


3.6955 


0.7177 


2s2p 


nS ^S^ 


0.0707 


0.7884 


2s2p ip° 


2p2 


9.2975 


10.0859 


2s2p ip° 


2s3d ^D^ 


3.7574 


13.8433 


2s2p ip° 


nZ) ^D^ 


2.4125 


16.2558 


2s2j3 ip° 


2p3p ip^ 


0.1613 


16.4171 


2s2p ip° 


nP ^P"" 


0.1169 


16.5340 


Core 




0.01964 


16.5536 



resonant transition and much of the difference with the CI+MBPT calculation comes from 
this transition with the CI+MBPT calculation giving 8.918 a|] [13]. The CI+MBPT calcu- 



lation overestimated the 25^ ^S^ - 2s2p transition energy difference by 0.27%. So it is 
possible that part of the discrepancy with the CI+MBPT calculation could be removed by 
using the experimental energy difference when calculating the polarizability of the resonant 
transition. It should be noted that experimental energy differences were used in a recent 
CI+MBPT calculation of the polarizability of the 3s^ state of Si^"*" js^] . 

The present CICP calculation of the ground state polarizability does not take into con- 
sideration the contribution from the 25^ ^S*^ — >■ 2s2p ^P° transition. The oscillator strength 
for this transition is only 3.361 x 10~^ jloj so this transition can be safely omitted from 
the determination of the polarizability. This also justifies the omission of the spin-orbit 
interaction from the effective Hamiltonian for the valence electrons. 
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D. The BBR shift 



The blackbody radiation shift of an atomic clock transition can be approximately calcu- 
lated using the Eqs. ([1]) and ([2]). In this expression the temperature in K is multiplied by 
3.1668153 X 10~^. Using the present polarizabilities and converting to frequency shifts at 
300 K gives Az/2,2 13.= = -0.08305 Hz and Az/2.2p spg = -0.06699 Hz. In the present CICP 
calculation the dipole polarizability difference for the 2s^ ^S*^ — )■ 2s2p ^Pq clock transition is 
Aai = -1.8643 al- The relativistic CI+MBPT calculation gave Acti = -1.851 ajj. 

Using a value of Aai = —1.8643 leads to a net frequency shift at 300 K of Ai/ = 0.01605 
Hz. This is consistent with the CI+MBPT result Az/ = 0.0159(16) Hz [isj]. A small 
correction to the polarizabilities needs to be made to allow for the slight variation of the 
polarizabilities due to the finite temperature of the BBR radiation field, 

ai(T) = ai(l + r/), (15) 

where ai{T) is the polarizability after correction, vj is the dynamic correction factor. The 
leading order term of t] is given by [12] 

407r2T2 

The value of 77 was found to be quite small. In the present CICP calculation, it was —1.42 x 
10~^ for the 2s^ ^S^ state and —1.45 x 10~^ for the 2s2p ^Pq state. Taking this correction 
into account, the 300 K dipole polarizabilities of 2s^ ^S^ state and 2s2p ^V" state are 9.6428 
Og and 7.7787 Cq respectively. The polarizability difference is Aai = —1.8642 ajj. This is 
only 0.0001 Oq smaller than the T = K value of —1.8643 Oq. It is evident that the effect of 
the dynamic correction in the B+ 2s^ ^S*^ — )■ 2s2p '^Pq clock transition is miniscule. 

When the uncertainties in the polarizabilities are taken into consideration the final recom- 
mended CICP polarizability difference at 300 K is —1.86(6) Og. The derived frequency shift 
is 0.0160(5) Hz. The CI+MBPT calculation gave a polarizability difference of —1.85(18) a\ 
and a frequency shift of 0.0159(16) Hz. The difference between the CICP and CI+MBPT 
calculations of the frequency shift is less than 1.0%. 

The uncertainty associated with the CI+MBPT calculation is more than three times 
larger than that quoted for the present CICP calculation. Although uncertainties are not 
assigned to the CI+MBPT polarizabilities, their final BBR shift uncertainty indicates un- 
certainties in their polarizabilities of 1.0%. The CI+MBPT uncertainty estimates seem very 
conservative given the 0.02 level of agreement between the CICP and CI+MBPT polar- 
izabilities. A more recent CI+MBPT calculation of the polarizability of the Si^"^ ground 
state quoted an uncertainty of 0.12% (sTj . 



IV. CONCLUSIONS 

The polarizabilities of some low lying states of the B^"*" and B+ ions are computed with 
large basis CI calculations with an underlying semi-empirical Hamiltonian. The motivation 
for these calculations was an independent calculation of the BBR shift of the B"*" 2s^ ^S*^ — )■ 
2s2p '^Vl clock transition [13]. 

The final estimate of the frequency shift, namely 0.0160(5) Hz is within 1% of the earlier 
CI+MBPT calculations 3]. The almost perfect agreement between these two completely 
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independent calculations gives increased confidence in the respective reliabilities of both 
calculations. One reason for the good agreement between both calculations is that both 
calculations give very accurate solutions of the Schrodinger equation with respect to their 
underlying Hamiltonian. Both the CICP and CI+MBPT approximate the aspects of the 
physics, and in particular the core-valence interaction. The CICP calculation uses a HF 
plus semi-empirical polarization potential to simulate core-valence correlation effects. The 
CI+MBPT calculation uses MBPT to incorporate the dynamical effects going beyond the 
HF interaction. Making these approximations simplifies the calculation sufficiently to allow a 
close to numerically exact solutions of the Schrodinger equation for the two valence electrons. 

In addition to the dipole polarizabilities, the present model computes the quadrupole, 
octupole and non-adiabatic dipole polarizabilities. One way to measure the B"*" polarizability 
would be the RESIS technique jsl, 59l|. The analysis of the raw experimental RESIS data 



can be improved if estimates of the quadrupole and non-adiabatic dipole polarizability are 
available. 
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